UNIVERSITY  OF  MARYLAND 
COMPUTER  SCIENCE  CENTER 

COLLEGE  PARK,  MARYLAND 
20742 

80  10  9  991 


t.  ! 


JR- 929  _ 

/aFCTsR-77-3271; 


LL > 


Aug«Rt , 


>  & 


iLl  CONVEX  DIGITAL  SOLIDS  #  '  (J 

■  J  ■?'  Chul  E .  '  Kiin  "\ 

Department  o I  Computer  Science 
University  of  Maryland 
College  Park,  MD  4G-T42 

Azriel yRosenfeld  ' 
Computer'  Vision  Laboratory 
Computer  Science  Center 
University  of  Maryland 
College  Park,  MD  20742 


.#  i 


ABSTRACT 

A  definition  of  convexity  of  digital  solids  is  introduced. 
Then  it  is  proved  that  a  digital  solid  is  convex  if  and  only 
if  it  has  the  chordal  triangle  property.  Other  geometric 
properties  which  characterize  convex  digital  regions  are 
shown  to  be  only  necessary,  but  not  sufficient,  conditions 
for  a  digital  solid  to  be  convex.  An  efficient  algorithm 
is  presented  that  determines  whether  or  not  a  digital  solid 
is  convex. 
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1. 


Introduction 


Discrete  geometry  lias  applications  in  image  processing, 
pattern  recognition  and  other  areas.  Before  an  image  of  an 
object  is  processed  by  a  computer,  it  is  digitized  to  yield  a 
finite  subset  of  digital  points  called  a  digital  image.  Then, 
various  types  of  operations  such  as  template  matching  and  geo¬ 
metric  property  measurement  are  performed  on  sets  of  digital 
points . 

Convexity  is  an  important  and  useful  geometric  property. 
Since  digital  image  processing  has  been  restricted  almost  en¬ 
tirely  to  2-dimensional  images,  digital  convexity  has  been 
studied  only  lur  digital  regions,  tnat  is,  8-connected  [10] 
finite  subsets  of  digital  points  in  the  plane.  Sklansky  in 
[13]  defined  a  digital  region  to  be  convex  if  and  only  if  there 
is  a  convex  region  wnose  image  (under  digitization)  is  the  digi¬ 
tal  region.  He  men  showed  that  a  digital  region  is  convex  if 
and  only  if  its  minimum-perimeter  polygon  is  convex  [13,14].  A 
digital  region  i<  is  said  to  have  the  line  property  if  every 
digital  point  on  the  line  segment  between  any  two  points  of  R 
ls  a  point  of  R.  Ln  17]  Minsky  and  Papert  defined  a  digital 
region  to  be  convex  if  and  only  if  it  has  the  line  property. 
However,  no  useful  results  regarding  convex  digital  regions  have 
been  derived  from  this  definition.  Let  d^,d^  be  two  points  of  R 
and  R'  the  union  of  all  cells  (grid  squares)  whose  centers  are 
points  ot  R.  Then  I’(R;d^,d^)  denotes  the  area  bounded  by  the 
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line  segment  d-j^  and  the  boundary  of  R1.  R  is  said  to  have 
the  area  property  if  for  any  two  points  d^,d£  of  R,  every  dig¬ 
ital  point  of  P(R;d^,d2)  is  a  point  of  R.  In  [2],  a  digital 
region  is  defined  to  be  convex  if  and  only  if  it  has  the  area 
property.  Moreover,  it  is  shown  in  [2,4]  that  the  above  three 
definitions  are  equivalent,  the  definition  of  Sklansky  having 
been  slightly  modified. 

In  [11]  Rosenfeld  introduced  the  chord  property  to  charac¬ 
terize  digital  straight  lines.  The  chord  property  and  the  area 
property  turned  out  to  be  equivalent  [5] .  Furthermore,  it  was 
shown  in  [3,4]  that  a  digital  region  is  convex  if  and  only  if 
any  two  points  of  the  region  may  be  connected  by  a  digital 
straight  line  segment  in  the  region. 

These  results  lead  us  to  believe  that  convexity  of  digital 
regions  is  now  well  defined  and  understood.  Because  of  the 
recent  growing  interest  in  3-dimensional  image  processing 
[1,6,9,16],  it  seems  essential  to  develop  a  theory  of  3-dimen¬ 
sional  discrete  geometry. 

In  this  paper  we  introduce  a  definition  of  convexity  of 
digital  solids,  which  is  an  extension  of  that  of  Sklansky.  Then 
we  show  that  a  digital  solid  is  convex  if  and  only  if  it  has 
the  cnordal  triangle  property.  However,  this  characterization 
of  convexity  does  not  lend  itself  to  development  of  an  efficient 
algorithm  to  recognize  convex  digital  solids.  Thus,  convexity 
of  digital  solids  is  further  characterized  in  terms  of  semi- 


digital  points  on  the  surface  of  the  convex  hull  of  the 
digital  solid.  With  this  characterization,  an  efficient 
algorithm  is  developed  to  determine  whether  or  not  a  given 
digital  solid  is  convex. 

The  geometric  properties  that  characterize  convex  digital 
regions  are  either  well  defined  in  or  easily  extendable  to  the 
3-dimensional  case.  It  will  be  shown  that  the  chordal  triangle 
property,  which  is  an  extension  of  the  chord  property,  is  the 
only  one  that  is  useful  for  characterizing  the  convexity  of 
digital  solids. 

In  the  next  section,  we  introduce  a  few  relevant  defini¬ 
tions,  notation  and  terminology.  Section  3  introduces  a  defi¬ 
nition  of  convexity  of  digital  solids  and  discusses  its  proper¬ 
ties.  An  efficient  algorithm  for  recognizing  convex  digital 
solids  is  presented  in  Section  4.  The  next  section  is  con¬ 
cerned  with  geometric  properties  that  characterize  convex  digi¬ 
tal  regions  and  solids. 


2. 


Definitions 


Consider  the  set  of  all  lattice  points  in  3-dimensional 
Euclidean  space.  A  lattic  point  d  =  (h,k,m)  is  called  a  digi¬ 
tal  point.  To  each  digital  point  d  is  associated  a  unit  cube 
whose  center  is  the  digital  point  and  faces  are  parallel  to 
the  coordinate  planes.  The  cube  associated  with  digital  point 
d  is  denoted  by  c  and  is  called  a  cell.  Digital  point  d'  is 
a  6-,  18- ,  or  26-neighbor  of  digital  point  d  if  c'  shares  with 
c  a  face,  an  edge  or  a  corner  point,  respectively. 

Let  R  be  a  set  of  digital  points.  Then  R  denotes  its  com¬ 
plement,  th 't  is,  the  set  of  all  digital  points  that  are  not 
in  R.  We  denote  by  R*  the  set  of  all  cells  that  are  associated 
with  the  points  in  R.  The  set  of  points  cf  3-dimensional  Eucli¬ 
dean  space  in  R1  is  denoted  by  s (R) ,  and  its  boundary  surface 
by  3s (R) .  For  any  two  points  x  and  y,  xy  denotes  the  line  seg¬ 
ment  between  the  two. 

A  chain  is  a  finite  sequence  of  digital  points  such  that 
every  element  of  the  sequence  except  the  first  is  a  6-neighbor 
of  its  predecessor.  A  set  R  of  digital  points  is  said  to  be 
6-connected  if  for  any  two  points  d,d'  of  R,  there  is  a  chain 
from  d  to  d'  in  R. 


Digital  solids 


A  digital  solid  S  is  a  finite  set  of  digital  points  which 


is  6-connected.  A  digital  solid  S  is  said  to  be  simple  if 


there  exists  no  pair  of  points  d^,d^  s  such  that  d^d^,  is 
parallel  to  a  coordinate  axis  and  there  is  a  point  of  S  on 
d^d^*  We  denote  H(S)  the  convex  hull  of  S,  that  is,  H(S) 
is  the  smallest  convex  polyhedron  that  contains  S.  Note  that 
the  vertices  of  H(S)  are  points  of  S. 

Digital  image 

A  digital  solid  S  is  the  digital  image  of  a  solid  q,  denoted 
S=I  (q) ,  if 

(i  )  qi-s  (S)  ,  and 

(ii)  if  c  is  an  element  of  S',  then  c°fq^0,  where  c°  is 
the  interior  of  c. 

Half-cell  expansion  [15] 

Let  S  be  a  digital  solid  and  Q  the  set  of  all  corner  points 
of  the  cells  of  S'.  Tnen  Q  is  a  digital  solid,  the  points  of 
Q  being  considered  as  digital  points,  and  we  denote  Q  by  E(S), 
the  half-cell  expansion  of  S. 

Di c^i tal  convexity 

A  digital  solid  S  is  said  to  be  digitally  convex  (or  simply 
"convex")  if  it  is  simple  and  there  is  a  convex  solid  q  such  that 
E(S)=I(q),  that  is,  the  half-cell  expansion  of  S  is  the  image 
of  the  convex  solid  q. 

Chord a 1  _t ria nglc  property 

Let  w=(x,y,z)  and  w'= (x' ,y ' , z ' )  be  two  points  in  3-dimensional 


Euclidean  space. 


The  two  points  are  said  to  be  near  each  other 


if  max{ | x-x' | , | y-y ' | | z-2 ' | } <1 .  Note  that  no  two  digital 
points  are  near  each  other. 


Let  S  be  a  digital  solid.  A  point  w  is  said  to  be  near 
S  if  there  is  a  point  of  S  which  w  is  near.  Let  d^,d2  and 
d^  be  points  of  S,  which  are  not  necessarily  distinct.  The 
triangle,  possibly  a  degenerate  one,  whose  vertices  are  d^, 
d2  and  d^,  is  called  a  chordal  triangle  of  S.  A  chordal  tri¬ 
angle  T  is  said  to  lie  near  S  if  every  point  of  T  is  near  S. 

A  digital  solid  S  is  said  to  have  the  chordal  triangle  property 
if  every  chordal  triangle  of  S  lies  near  it. 


The  main  result  of  this  section  is  that  the  chordal  triangle 
property  is  a  necessary  and  sufficient  condition  for  a  digital 
solid  to  be  convex.  First  we  show  that  a  digital  solid  S  has 
the  chordal  triangle  property  if  and  only  if  every  point  of 
H(S),  the  convex  hull  of  S,  is  near  S.  Next  it  is  shown  that 
every  point  of  H(S)  is  near  S  if  and  only  if  S  is  digitally 
convex.  The  main  result  then  follows  immediately. 

Lemma  1 :  A  digital  solid  S  has  the  chordal  triangle  property 
if  and  only  if  every  point  of  H(S)  is  near  3. 

Proof :  Suppose  tnat  every  point  of  H(S)  is  near  S,  and  let  T 

be  a  chordal  triangle  of  S.  Since  T  is  a  subset  of  H(S),  every 
point  of  T  is  near  S  and  T  lies  near  S.  Therefore,  S  has  the 
chordal  triangle  property. 

Now  suppose  that  there  are  points  of  H(S)  which  are  not 
near  S,  ana  let  w  be  such  a  point.  If  w  is  on  a  face  of  H(S) , 
then  obviously  w  is  a  point  of  a  chordal  triangle.  Hence,  S 
does  not  have  the  chordal  triangle  property.  Assume  that  w  is 
an  interior  point  of  H(S).  We  further  assume  that  w  is  not  a 
digital  point.  (It  turns  out  that  the  case  where  w  is  a  digital 
point  is  taken  care  of  when  we  consider  the  case  where  w  is  not 
a  digital  point.)  Let  c  be  the  cell  of  which  w  is  a  point.  Then 
the  digital  point  d  which  is  the  center  of  c  is  not  a  point  of 
S,  since  otherwise  w  is  near  S.  There  are  two  cases  to  consider 


Case  1 : 


The  digital  point  d  is  not  a  point  of  H(S). 


Let  u  be  the  point  at  which  dw  intersects  a  face  of  H(S). 

We  claim  that  u  also  is  not  near  S.  If  so,  u  is  a  point  on  a 
face  of  H(S)  and  not  near  S.  Thus,  S  does  not  have  the  chordal 
triangle  property.  It  remains  to  prove  our  claim.  If  d=(h,k,m) 
then  w=(h+Ax,  k+Ay,  m+Az),  where  |Ax|*l/2,  | Ay | *1/2  and  |  Az  |  *1/2 

Without  loss  of  generality  assume  that  0 *Ax , Ay , Az *1/2 .  Then 
the  point  u  is  such  that  u=  (h+Ax1  ,k+Ay ' ,m+Az 1 )  and  0*Ax'*Ax, 
0*Ay'*Ay  and  0*Az'*Az.  (See  Figure  1.)  Suppose  u  is  near  S. 


Then  there  is  a  point  d.  in  S  which  is  near  u.  If  d,  =  (h+<$  ,k+5 

1  1  x  y 

m+^z),  tnen  S  =1  if  At>0  and  (S t~0  if  At=0  for  all  t=x,  y  and  z. 
Thus  w  is  also  near  d^,  which  is  a  contradiction.  This  proves 
our  claim. 

Case  2:  The  digitai  point  d  is  a  point  of  H(S). 

Obviously,  d  is  not  near  S  (if  it  were,  it  would  have  to  be 
in  S,  making  w  near  S) .  As  mentioned  above,  this  also  takes 
care  of  the  case  where  w  is  a  digital  point  (so  that  w=d) . 


Consider  the  set  X  of  all  digital  points  of  H(S)  that  lie  on 
the  line  t  which  passes  through  d  and  is  parallel  to  the 
x-axis.  If  X  contains  points  of  S  on  both  sides  of  d,  let  u 
and  v  be  the  nearest  points  of  S  on  each  side  of  d.  The  line 
segment  uv  is  a  degenerate  chordal  triangle  of  S  and  it  does 
not  lie  near  S.  Therefore,  S  does  not  have  the  chordal  triangle 
property.  Suppose  that  X  does  not  contain  any  point  of  S  on 
one  side  of  d.  Let  u  be  the  point  at  which  l  intersects  a 
face  of  H(S)  and  there  is  no  point  of  S  on  du.  We  claim  that 
u  is  not  near  S  and  thus  S  does  not  have  the  chordal  triangle 
property.  Now  we  prove  our  claim.  If  u  is  a  digital  point, 
then  obviously  it  is  not  near  S.  So  assume  that  u  is  not  a 
digital  point.  If  d=(h,k,m),  then  u=(n+x,k,m)  because  u  is 
a  point  on  Also  x  is  not  an  integer,  since  u  is  not  a  digi¬ 

tal  point.  The  only  digital  points  near  u  are  d ’ =  (h+ t x j  , k , m) 
and  d"= (h+ f x 1 , k , m) ,  where  lx]  is  the  largest  integer  not 
greater  than  x  and  fxl  the  smallest  integer  not  less  than  x. 
Neither  d'  nor  d"  is  a  point  of  S  and  u  is  not  near  S.  So 
our  claim  is  proved.  This  completes  the  proof  of  this 
lemma .  L. 

We  need  a  few  lemmas  before  we  derive  a  result  on  the  rela¬ 
tionship  between  the  convex  hull  and  convexity  of  a  digital  solid. 
Lemma  2 :  Let  S  be  a  digital  solid.  If  c  is  a  cell  of  E(S) , 
then  cLH  (S)  /0 . 

Proof :  If  c  is  a  cell  of  E(S),  at  least  one  corner  point  of  c, 
say  e,  is  a  point  of  S.  Then  e eel  ill  (S).  l.. 


Lemma  3 :  Let  S  be  a  digital  solid.  If  every  point  of  H  (S)  is 
near  S,  then  S  is  simple. 

Proof :  Suppose  every  point  of  H(S)  is  near  S.  Then  every 
digital  point  of  H (S)  is  a  point  of  S.  Consider  any  pair  of 
points  d^,d2  of  S  such  that  is  parallel  to  a  coordinate 

axis.  Since  H(S)  is  convex,  d^d2  is  a  subset  of  H(S)  .  There¬ 
fore,  every  diyital  point  on  <3^2  is  a  point  of  S  and  S  is 

simple  by  definition.  L 

Lemma  4 :  Let  S  be  a  digital  solid.  If  every  point  at  H(S)  is 
near  S,  then  H (S)  is  a  subset  of  the  set  of  interior  points  of 
s  ( L  ( S ) ) .  Therefore,  there  exists  a  positive  number  s  such  that 
dist(u,v)*e  for  any  point  u  of  H(S)  and  any  point  v  of  ds(E(S)). 
Proof :  Let  w  be  a  point  of  3s(E(S))  and  c  be  a  cell  of  E(S) 

to  which  w  belongs.  If  w  is  a  corner  point  of  c,  then  it  is 
not  near  S.  If  w  is  on  an  edge  of  c,  neither  end  point  of  the 
edge  is  a  point  of  S ;  hence  w  is  not  near  S.  If  w  is  on  a  face 

of  c,  no  corner  point  of  the  face  is  a  point  of  S,  and  again  w 

is  not  near  S.  Therefore,  H(S)  and  3s(E(S))  do  not  meet  or 
intersect.  Since  every  point  of  S  is  an  interior  point  of  s(E(S)) 
the  vertices  of  H(S)  are  interior  points  of  s(E(S)).  Thus,  H(S) 
is  a  subset  of  the  set  of  interior  points  of  s(E(S)). 

Since  dS(E(S))  and  H(S)  are  both  compact  sets  and  do  not 

have  any  point  in  common,  min  .  ..  dist  (u ,  v)  =e  >0  .  [_. 

■*  ufc.ll  (S)  ,  vfcds  (E  (S)  ) 


*  tr- 


Lemma  5 :  Let  S  be  a  digital  solid.  If  q  is  a  convex  solid  such 
that  E(S)sl(q),  then  H(S)^q°,  where  q°  is  the  interior  of  q. 
Proof :  Let  v  be  a  vertex  of  H(S) .  Then  there  are  eight  cells 
of  E(S)  of  which  v  is  a  corner  point.  Since  H(S)  is  convex, 
there  is  one  cell  among  the  eight,  say  c,  such  that  c°riH  (S )  =0 . 
Therefore,  s  (E  (S)  )  (H  (S)  )  .  Since  obviously  H(S)s-q,  H(S)  must 
be  a  proper  subset  of  q.  Suppose  r=3qli8H(S),  the  intersection 
of  the  boundary  surfaces  of  q  and  H(S),  is  not  empty.  Since  q 
is  convex,  r  contains  a  vertex  of  H(S),  say  v,  and  there  is  a 
plane  p  such  that  v  is  on  p  and  q  lies  in  one  side  of  p.  Then 
there  is  a  cell  of  E(S),  call  it  c,  such  that  v  is  a  cornet- 
point  of  c  and  c  lies  on  the  opposite  side  of  p  from  q.  Thus, 
cofq=0  and  E(S)^I(q),  which  is  a  contradiction.  Therefore,  r 
is  empty  and  H  (S )  '-q° .  L 

Lemma  6 :  A  digital  solid  S  is  convex  if  and  only  if  every 
point  of  H (S )  is  near  S. 

Proof :  Suppose  that  every  point  of  L(S)  is  near  S.  Let  q  be 

a  convex  polyhedron  obtained  from  11  (S)  by  parallel  translation 
of  all  its  faces  outward  by  a  distance  .  By  Lemma  4,  there 
exists  a  positive  number  ,  which  depends  on  t  in  the  lemma, 
such  that  q  is  a  subset  of  s(E(S)).  If  c  is  a  cell  of  E(S), 
chH(S)^0  by  Lemma  2  and  thus  c°iq^0  since  q  is  an  expansion  of 
H  ( S ) .  Hence,  E(S)=I(q).  Also  by  Lemma  3,  S  is  simple.  There¬ 
for,  S  is  digitally  convex. 


Now  suppose  that  there  are  points  of  H(S)  which  are  not 


near  S.  Since  every  interior  point  of  s(E(S))  is  near  S, 

H  (S )  contains  a  point  w  which  is  either  a  point  of  9s(E(S))  or  of 
s(E(S)).  Thus,  w  is  a  point  of  a  cell  c  which  is  not  a  cell  of 
E (S ) .  Suppose  there  exists  a  convex  solid  q  whose  image  is 
E(S).  Since  E(S)s-l(q),  H(S)£~q°  by  Lemma  5.  Thus  w  must  be  an 
interior  point  of  q  and  c°fq^0.  Then  ctl (q)  but  c  is  not  a 
cell  of  E(S),  which  is  a  contradiction.  Therefore,  S  is  not 
digitally  convex  because  there  is  no  convex  solid  whose  image 
is  E  (S)  .  L 

The  main  result  of  this  section  follows  immediately  from 
Lemmas  1  and  6: 

Theorem  7 :  A  digital  solid  is  convex  if  and  only  if  it  has 
the  chordal  triangle  property. 


4. 


In  the  previous  section  we  showed  that  the  chordal  triangle 
property  is  a  simple  geometric  property  that  characterizes 
convex  digital  solids.  However,  it  does  not  lend  itself  to 
development  of  an  efficient  algorithm  to  determine  whether  or 
not  a  digital  solid  is  convex.  The  reason  is  that  there  are 
infinitely  many  points  on  a  chordal  triangle  and  the  number  of 

chordal  triangles  is  O(n^),  where  n  is  the  number  of  points  in 

S.  In  the  sequel,  it  will  be  shown  that  only  finitely  many 
points  on  the  surface  of  H(S)  need  be  examined  to  recognize 
the  convexity  of  a  digital  solid. 

Let  p  be  a  polygon  on  a  plane  in  3-dimensional  Euclidean 
space.  A  point  w=(x,y,z)  of  p  is  said  to  be  a  semi-digital 
point  if  at  least  two  of  its  coordinates  are  integers  or  it  is 
on  an  edge  of  p  and  one  of  its  coordinates  is  an  integer.  Note 
that  the  number  of  semi-digital  points  on  p  is  approximately 
equal  to  its  area. 

Our  algorithm  is  based  on  the  result  stated  in  the  following 
theorem. 

Theorem  3:  A  simple  digital  solid  S  is  convex  if  and  only  if 

every  semi-digital  point  on  the  faces  of  H(S)  is  near  S. 

Proof :  "Only  if"  is  a  special  case  (if  Lemma  6.  Conversely, 
suppose  that  the  simple  digital  solid  S  is  not  convex.  By 
Lemma  0,  there  arc  points  of  H(S)  which  are  not  near  S;  let 


w=(x,y,z)  be  such  a  point.  There  are  two  cases  to  consider. 

Case  1 :  w=(x,y,z)  is  a  digital  point,  i.e.,  x,y,z  are  all 

integers . 

Let  t  be  the  ray  from  w  parallel  to  the  x-axis;  then  i 
contains  no  point  of  S  because  S  is  simple.  Let  u  be  the 
point  at  which  Z  intersects  a  face  of  H(S).  Then  u  is  not 
near  S  and  u=(x',y,z)  is  a  semi-digital  point. 

Case  2:  w=(x,y,z)  is  not  a  digital  point. 

Let  c  be  the  cell  of  which  w  is  a  point.  Then  d,  the 
center  of  c,  is  not  a  point  of  S,  since  otherwise  w  is  near  S. 

[f  d  is  a  point  of  H(S),  then  Case  1  applies  to  it,  and  there 
is  a  semi-digital  point  on  a  face  of  H(S).  Thus,  assume  that  d 
is  not  a  point  of  H(S).  Let  u  be  the  point  at  which  dw  inter¬ 
sects  a  face  of  H(S)  (see  Figure  1.).  Then  u  is  not  near  S. 
(ReLer  to  the  proof  of  Lemma  1  for  the  reason.)  Let  d=(h,k,m). 
Then  u=  (U  *Ax ,  k.+Ay ,  m+Az )  ,  where  |  Ax  |  **1/2  ,  |  Ay  |  **1/2  and  |Az|*»l/2. 

(a)  u  is  on  an  edge  of  H  (S)  . 

As  we  move  u  along  the  edge  in  either  direction.  Ax, Ay  and 
Az  change  their  values  continuously  and  monotonically .  Let  us 
move  u  until  one  of  Ax, Ay  and  Az  attains  a  value  of  0,  1,  or  -1 
for  tlie  first  time.  This  always  occurs  because  the  endpoints  of 
the  edge  are  digital  points.  The  resulting  u  is  a  semi-digital 
point.  It  is  easy  to  see  that  u  is  not  near  S. 

(b)  u  is  not  on  an  edge  of  H(S). 

Take  a  coordinate  plane  which  is  not  parallel  to  the  face, 
say  the  xy-plane.  Consider  the  line  t  of  intersection 

Xy 


xy 


of  Liu;  lace  and  the  plane  which  is  parallel  to  tne  xy-pi  am- 
and  contains  the  point  u.  As  the  point  u  is  moved  along  i 
tne  values  of  Ax  and  .Ay  change  wniLe  that  of  ./  is  Kept  con¬ 
stant.  We  move  u  until  either  it  reaches  an  e<i  ol  the  face 

or  one  of  Ax  and  Ay  attains  a  value  o!  0,  1,  m  -]  tor  the 

first  time.  It  u  reaches  an  edge  first,  tiien  we  have  Case  Jta)  , 
which  has  already  been  treated.  Now  assume  without  loss  ot 
generality  that  Ax  attains  a  value  of  0,  1,  or  -i  first.  Thou 

the  x-eoot'di  nate  of  u  is  an  integer.  Cons  ulo:  tne  i  i  n«  •  t  <  ■ : 
intersection  of  tne  face  and  the  plane  which  is  out ai  lei  f.  t ne 
yz-plane  and  contains  u.  (If  tne  //.-piano  wet  o  para  I  It;  1  t  •  >  *  no 
fact',  then  ’  y  would  have  attained  a  value  ol  0,  1,  m  i  :  r  . :  *  , 

since  Ax  would  have  stayed  constant.  In  this  case  i  w  e;id  ne 

x  z 

considered  instead.)  Move  tne  point  u  along  t  ,  until  eithet  it 

reaches  an  edge  ot  the  face  or  the  value  ot  y  or  '/  Decomes  0  , 

1,  or  -1  lor  the  first  time.  If  u  reaches  an  odqe  first,  ttien 
it  is  a  semi-diyitai  point.  If  Ay  or  Az  attains  a  value  of  0, 

1  or  -1  first,  then  again  u  is  a  semi-diyitai  point,  since  two 
of  its  coordinates  are  inteyors.  It  can  easily  be  shown  that  u 
is  not  near  b.  L. 

Corollary  9;  A  digital  solid  S  is  convex  if  and  only  if  it  is 
simple  and  every  semi-digital  point  on  the  surface  of  H(S)  is 
near  S. 

We  are  now  ready  to  present  an  algorithm  to  determine  whethei 


or  not  a  given  digital  solid  is  convex.  In  fact,  Corollary  9 


concisely  describes  the  algorithm.  The  algorithm  first  checks 
il  S  is  or  is  not  simple.  If  S  is  not  simple,  then  the  algo¬ 
rithm  outputs  that  S  is  not  convex  and  halts.  Otherwise  it 
examines  every  semi-digital  point  on  the  surface  of  H(S)  to 
see  whetnei  any  one  ot  them  is  not  near  S.  If  such  a  point  is 
found,  trie  algorithm  halts  after  it  outputs  that  S  is  not  con¬ 
vex.  it  all  semi-digital  points  are  near  S,  then  S  is  convex. 
Below  t ue  ulgontnm  is  presented  formally. 

Algol  1  thm  JO-CONVI.XITY  (S) 

I.  If  S  is  not  simple  then  output  (False ) ;  stop. 

J.  Construct  the  convex  hull  H ( S )  of  S. 


1.  It  a  serai-digital  point  which  is  not  near  S  is  found 
on  the  surface  of  H(S)  then  output  (False) ;  stop. 

4.  output  (True);  stop. 

The  correctness  ot  the  algoritnm  is  immediate  from  Corollary 
).  However,  determination  of  its  computational  complexity  re¬ 
quires  a  detailed  description  of  each  step  and  the  data  struc¬ 
ture:.  used  in  tlie  algorithm.  For  simplicity  we  assume  that  a 
digital  solid  S  is  a  subset  of  the  set  of  n^  digital  points  in 
tne  cube  whose  edge  is  ot  length  n.  S  is  represented  by  a  run 
length  code  112|  such  that  RC(i,j)  is  a  finite  sequence  of  run 
legnths  of  U's,  digital  points  of  S  with  coordinates  (i,j,z), 
and  l's,  digital  points  of  S  with  coordinates  (i,j,z).  Thus, 


RC  ( i  ,  i )  (f  ^  y  ,  f  x  ■  j,  .  .  .  ,  l  .  .  )  represents  the  (i  ,  j  )  th  row  as 

M  ljrij 

composed  of  a  run  of  0's  of  length  followed  by  a  run  of  l's 

of  length  t .  . .  and  so  on. 


and  so  on. 


.✓  • 


( I )  Is  S  simple? 

11  a  row  has  more  than  one  run  of  l's,  then  S  is  obviously 
not  simple.  Consider  the  set  of  rows  in  a  piano  parallel  to  the 


yz-plane,  that  is,  the  set  of  rows  (i,j)  for  all  3,  1  a»j  s-n ,  and 

a  fixed  i,  lM^n.  Suppose  that  f.  ...  M.  .  ,  A  and  t’  n>f.  ,  ,  .  n 

^  13O  i,]  +  l, 0  ikO  i,k  +  l,0 

for  some  lij-sk^-n.  Then  S  is  not  simple.  Also,  if  £ .  .„+!’.  .  > 

J  r  l  j  0  l  j  1 

(  .  .  ,  ,  ..  +  !'.  ■  ,  ,  .  and  1’  .  ,  ,  „  +  i  .  ,  ,  ,  ,  ,  then  S  is  not 

l, 3+1,0  i, 3+1,1  ikO  lkl  i,k+l,0  i,k+l,l 

simple.  The  set  of  rows  in  a  plane  parallel  to  the  xz-plane  can 
be  checked  similarly.  It  is  easy  to  see  that  S  is  simple  other¬ 
wise.  The  above  observations  lead  to  algorithm  SIMPLE  for  step 


1  of  algorithm  JU-CONVEX1TY . 


Algorithm  SIMPLE (S) 

1.1.  For  each  i  and  eacn  j,  ivi,j*n, 

check  RC(i,j)  to  see  if  there  is  more  than  one 


run  of  l's;  if  so,  output  (False) ;  stop. 

1.2.  For  each  i,  l*i*n, 


check  if  1  „  increases  and  then  decreases  or 

13O 

<  .  „+£  . ,  decreases  and  then  increases  as  3  increases 

13O  13I  J 

from  1  to  n;  if  so,  output  (False) ;  stop. 

1.3.  For  each  3,  l*j*n, 

check  if  l . . „  increases  and  then  decreases  or 
13O 

£ .  . _+£. decreases  and  then  increases  as  i  increases 
13O  13I 

from  1  to  n;  if  so,  output  (False)  ;  stop. 

2 

Each  step  runs  in  0(n  )  time  and  requires  a  constant  work 


space . 


(2)  Construction  of  H(S) 

A  point  of  S  is  a  corner  point  if  three  of  its  6-neighbors 

are  points  of  S  and  they  are  mutually  18-neighbors.  Only  a  corner 

point  of  S  may  be  a  vertex  of  H(S),  and  the  convex  hull  of  the 

set  of  corner  points  of  S  is  also  the  convex  hull  of  S.  Obviously, 

only  the  first  and  the  last  digital  points  of  each  row  can  be 

2 

corner  points.  Hence,  there  are  at  most  0(n  )  corner  points 

2 

in  S  because  there  are  n  rows.  Let  CP(S)  denote  the  set  of 
corner  points  of  S. 

An  algorithm  to  build  the  convex  hull  of  a  set  of  points  is 
given  in  [8].  We  denote  the  algorithm  by  HULL.  Its  input  is  a 
set  of  points  K  and  its  output  is  the  convex  hull  H (R) = (F^ , . . . , ) 
where  each  face  I'V  of  H  (R)  is  represented  by  a  sequence  of  its 
vertices . 

Algorithm  CONVLXHULL (S, H (S) ,  k) 

2.1.  Obtain  CP(S),  the  set  of  corner  points  of  S,  by  checking 
the  tirst  and  last  points  of  S  in  each  row  (i,j),  where 

1  *i  ,  ]  a-n  . 

2.2.  Call  HULL (CP(S) ,H(S) )  to  construct  H (S ) = (F^ , . . . , F^ ) . 

2.3.  Return . 

2 

Step  2.1  requires  0(n  )  computing  time  and  work  space.  Step 
2  2 

2.2  runs  in  0{n  log  n)  time  and  needs  0(n  )  work  space  to  store 


I  i )  Is  every  semi-diqital  point  ne.ir  S? 

Mi  nee  the  number  of  semi-diqital  points  is  approximately 

equal  to  tile  area  of  the  surface  of  11  (3)  and  H  ( S )  is  convex, 

> 

there  are  0(n  )  semi-diqital  points.  Given  a  face  of  H(S) 
represented  by  a  sequence  of  its  vertices,  each  semi-diqital 
point  may  be  located  in  constant  con, [Hit  in>:  time.  If  a  semi- 
diqital  point  is  a  diqital  [joint,  then  it  is  near  S  if  and  only 
it  it  is  a  point  of  S.  If  a  semi-diqital  point  has  two  inteqer 
coordinates,  then  two  Jiqitai  points  must  be  checked.  For  in¬ 
stance,  if  a  semi-diqital  point  is  w-  (ii  +  .'.x ,  k  ,  m)  ,  0  v'.x<l ,  then 
it  is  near  S  it  either  (h,k,m)  or  (h+l,k,m)  is  a  point  of  S. 
it  a  semi-diqital  point  w  has  only  one  inteqer  coordinate,  then 
lour  diqital  [joints  must  be  checked  to  see  if  w  is  near  S. 

Thus,  we  have  tile  tollowinq  alqoi  l  turn  for  stup  3  of  alqorithm 
3!>-C<  iNVF.X  I  TY  . 

Algorithm  M.AR  (M ,  II  ( S )  ,  k ) 

3.1.  For  each  face  F  of  h(S),  l*i*k, 

for  each  seta -d l q i La  1  point  w  on  F 

•heck  if  w  is  iieai  S;  if  not,  out['ut  (False); 
s  to[  * . 

3.2.  Return . 

2 

obviously  alqoritnm  NFAR  has  tunnitnj  time  of  0(n  )  and  needs 


constant,  woik  space. 


Theorem  10:  Algorithm  3D-CONVEXITY  determines  whether  or  not 

2 

a  digital  solid  is  convex,  runs  m  0 (n  log  n)  time  and  requires 
2 

0(n  )  work  space. 


5. 


Convex  digital  regions  and  solids 


We  restate  formal  definitions  of  the  geometric  properties 
that  are  used  to  define  convexity  of  digital  regions. 

A  digital  region  R  is  said  to  have  the  line  property  if  there 
is  no  triplet  (d^d^d^)  of  collinear  digital  points  such  that 
d^  and  d^  are  points  of  R  and  d^  is  a  point  of  R. 

Area  property 

Let  ls(R)  denote  the  boundary  of  the  set  of  points  of  cells 
whose  centers  are  the  points  of  R.  R  is  said  to  have  the  area 
property  if  there  are  no  two  points  d^,d2  of  R  such  that  the 
bounded  area  whose  boundaries  consist  of  nonempty  segments  of 
d-^d'2  and  ds(R)  contains  a  digital  point  of  R. 

Chor d  property 

A  digital  region  R  is  said  to  have  the  chord  property  if 
for  each  point  w=(x,y)  on  d , d „ ,  where  d,  and  d0  are  any  two 
points  of  R,  there  is  a  point  d=(h,k)  of  R  such  that 
maxi  i  x-h  |  ,  |  y-k  |  !  *-l . 

All  three  properties  above  are  equivalent  in  that  a  digital 
region  is  convex  if  and  only  if  if  has  any  one  of  them.  There 
are  Euclidean  geometric  properties  of  which  the  line  property 
and  the  area  property  may  be  considered  digital  equivalents, 
but  this  is  not  the  case  for  the  chord  property.  It  is  obvious 
that  the  line  and  chord  properties  are  well  defined  in  3-dimen- 
sionaJ  discrete  geometry.  It  is  immediate  that  the  chordal 


triangle  property  is  an  extension  of  the  chord  property,  but 
there  is  no  natural  extension  of  the  line  property.  In  [2] 
it  was  shown  that  a  digital  region  R  has  the  area  property 
if  and  only  if  the  convex  hull  of  R  contains  no  point  of  Rj 
Thus,  the  volume  property  which  is  defined  below  may  be  con¬ 
sidered  as  an  extension  of  the  area  property. 

Volume  property 

A  digital  solid  S  is  said  to  have  the  volume  property  if 
H (S )  contains  no  point  of  S. 

In  Section  3,  we  proved  that  a  digital  solid  is  convex  if 
and  only  if  it  has  the  chordal  triangle  property.  In  the  sequel, 
we  show  that  every  other  property  mentioned  here  is  a  necessary 
but  not  a  sufficient  condition  for  a  digital  solid  to  be  convex. 
Theorem  11:  Each  one  of  the  line,  chord,  and  volume  properties 
is  a  necessary  but  not  a  sufficient  condition  for  a  digital  solid 
to  be  convex . 

Proof ;  Suppose  that  a  digital  solid  S  does  not  have  either  the 
line  property  or  the  chord  property,  then  trivially  S  does  not 
have  the  chordal  triangle  property.  Hence,  S  is  not  convex  by 
Theorem  7.  Now  suppose  that  S  does  not  satisfy  the  volume  pro¬ 
perty.  Then  not  every  point  of  H(S)  is  near  S  and  S  is  not  con¬ 
vex  by  Lemma  6. 


To  see  that  these  conditions  are  not  sufficient,  consider 
the  polyhedron  q  which  is  formed  by  the  three  coordinate  planes 


and  two  planes  whose  equations  are  x/3+y/3-l  and  x/3ty/6+z/2=l , 
respectively.  Let  be  the  digital  solid  that  consists  of  all 
digital  points  of  q  except  the  point  d= (1,1,1).  (Sec  Figure  2.) 


Figure  2.  A  digital  solid  which  is  not  convex  but  has 

both  the  line  property  and  the  chord  property. 

Then  q  is  li(S^),  the  convex  hull  of  S .  .  Since  it  is  not  a  point 

of  S  but  is  a  digital  point,  d  is  a  point  of  H(S^)  which  is  not. 

near  S^.  Thus,  is  not  convex.  Note  that  every  point  of  H(S^) 

except  d  is  near  S-^.  It  is  easy  to  see  that  there  are  no  two 

points  d^,d2  of  S  such  that  u  is  on  d^d^ .  Therefore,  has  both 

the  line  property  and  the  chord  property. 

Next  let  S2  be  the  digital  solid  { (0 ,0 ,0) , (0 ,0 , 1) , (0 , 0 , 2) , 

(0,1,2),  (1,1,2),  (1,0,0)  i.  (See  Figure  3.)  Then  the  point  w= 

(1, 1/2,1)  is  a  point  of  H  ( S 2 )  and  not  near  S2  ,  so  is  not  convex. 

Since  H(S2)  contains  no  point  of  S2 ,  S2  has  the  volume  property.  i__ 


TtWii 


Co,  I,  z) 


Figure  3.  A  digital  solid  which  is  not  convex  but  has 
the  volume  property. 


Conclusion 


(> . 

A  definition  of  convexity  of  digital  solids  was  introduced. 

It  is  a  3-dimensional  extension  of  the  modified  definition  of 
Sklansky  [3,13].  Just  as  convexity  of  digital  regions  is  charac¬ 
terized  by  the  chord  property,  convexity  of  digital  solids  is 
characterized  by  the  chordal  triangle  property,  which  is  a  3- 
dimensional  extension  of  the  chord  property.  The  main  question 
concerning  convexity  of  digital  solids  is  under  what  conditions 
a  digital  solid  is  the  digitization  of  a  convex  solid.  The 
definition  presented  in  this  paper  and  the  characterization  of 
convexity  by  the  chordal  triangle  property  seem  to  provide  a 
satisfactory  answer  to  the  question. 

There  are  other  geometric  properties  that  are  used  zo 
characterize  convexity  of  digital  regions.  Somewhat  surprisingly, 
they  or  their  3-dimensional  extensions  turn  out  to  be  only  neces¬ 
sary  conditions.  It  would  be  interesting  to  determine  the 
class  of  digital  solids  that  each  of  these  properties  charac¬ 
terizes  . 

If,  given  a  digital  solid,  a  sequential  algorithm  must  con¬ 
struct  its  convex  hull  to  determine  its  convexity,  the  time  com¬ 
plexity  obtained  here  cannot  be  improved.  For,  0 (k  log  k)  is  the 
optimal  time  complexity  for  any  sequential  algorithm  to  construct 
the  convex  hull  of  k  points.  Thus,  to  develop  a  faster  algorithm, 
convexity  of  digital  solids  must  be  characterized  by  a  simpler 
geometric  property. 
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